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We demonstrate how dynamic correlations of heavy-hole and light-hole excitons 
in semiconductor quantum wells may be investigated by two dimensional correlation 
spectroscopy (2DCS). The coherent response to three femtosecond optical pulses is 
predicted to yield cross (off-diagonal) peaks that contain direct signatures of many- 
body two-exciton correlations. Signals generated at various phase-matching direc- 
tions are compared. 

I. INTRODUCTION 

Understanding the signatures of many-body interactions in the nonlinear optical response 
of semiconductors is an important fundamental problem with implications to all-optical and 
electro-optical device applications. 1 The linear response to a weak optical field is well 
described by a model of non-interacting quasiparticles. However, residual interactions, 
not accounted for by these quasiparticles, can considerably affect the nonlinear response. 
Similar to Frenkel excitons in molecular crystals and aggregates^ 2 Coulomb correlations 
among quasiparticles can dominate the nonlinear optical response of semiconductors, in 
marked contrast to the behavior of atomic systems.-^^^&^ii 1 ^ 1 ^ 

The coherent ultrafast response and many-body correlations in semi- 
conductor heterostructures have been studied extensively in the past two 
decades.-^^^^^^^ 1 ^ 1 ^ 12 ^ 12 ^ 12 ^ Due to various dephasing and relaxation mecha- 
nisms, the coherent response usually persists only on the tens of picosecond time scale. 
Optical spectra such as the linear absorption, pump-probe and Four-wave mixing (FWM) 
are commonly displayed as a function of a single (time or frequency) variable, and hence 
provide a one-dimensional (ID) projection of the microscopic information. ID spectra 
are hard to interpret in systems with many congested energy levels. The spectroscopic 
signatures of complex many-body dynamics projected on a ID spectral plot strongly 
overlap and may not be easily identified. For example, when ID techniques are employed 
in III-V semiconductor quantum wells (QWs), it is difficult to pinpoint the signatures of 
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two-excitons due to the small (1 to 2 meV) two-exciton interaction energies, i.e., bound 
two-exciton binding energy and unbound two-exciton scattering energy.— i2L2& Two-exciton 
effects can be easily masked by homogeneous and inhomogeneous line broadening. When 
both light-hole (LH) and heavy-hole (HH) excitons and their interactions are taken into 
account, the situation becomes even more complicated since there are different types of 
two-excitons such as LH, HH and mixed two-excitons.— 1 ^ Even in II- VI semiconductor 
QWs, where the two-exciton interaction energies are much larger, ( from several to several 
tens of meV), 29 it is still hard to resolve them. 

Multi-dimensional spectroscopy^ 1 ^ 1 ^ 1 ^^ 4 - 1 ^ in which the optical response is recorded 
and displayed versus several arguments, can overcome these limitations by separating the 
signatures of different pathways of the density matrix, known as Liouville space pathways.— 
Because different pathways are usually connected with specific couplings, one can resolve 
different coherences and many-body interactions by focusing on different peaks in multi- 
dimensional correlation plots. Optical and infrared 2D Correlation Spectroscopy (2DCS) is 
a femtosecond analogue of mult i- dimensional nuclear magnetic resonance (NMR)^ 1 ^ that 
has been shown to be very powerful for probing the structure and dynamics in chemical and 
biological systems.— dMLMiMd&dliM These techniques can reveal the coupling strength of 
elementary excitation through cross-peaks. After several attempts 4 ^ 1 ^ 4 ^^ to go beyond 
ID techniques, the first experimental implementation of 2DCS to investigate the many-body 
Coulomb interactions among LH and heavy-hole HH excitons in semiconductors has been 
reported recently.— 1 ^ 

In this work, we present the basic principles^ 4 of 2DCS in semiconductor QWs and show 
how it can be used to study the couplings among LH and HH excitons in a GaAs single 
quantum well. The same methodology previously applied to Frenkel excitons (hard-core 
bosons) in molecular aggregates or anharmonic vibrations (soft-core bosons) in proteins^ 4 - 
can be extended to Wannier excitons.— <22i5ii5Si In Section II, we survey the possible third- 
order 2D techniques of excitons. The double-sided Feynman diagrams and the corresponding 
sum-over-states expressions for 2DCS in III-V semiconductor quantum wells are presented in 
Section III. These guide the analysis of the numerical calculations presented in later sections. 
In Section IV, we present the multiband Hamiltonian and the nonlinear exciton equations 
used for the numerical calculation of the 2DCS signals. The numerical calculations are 
analyzed in Section V. Finally, in Section VI, we summarize our findings. 
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II. LIOUVILLE-SPACE PATHWAY ANALYSIS FOR 2D CORRELATION 

SPECTROSCOPY OF EXCITONS 

Impulsive 2DCS signals generated by well-separated femtosecond pulses may be calculated 
using sum-over- states expressions which provide insights into the origin of various peaks 
associated with different coherences and many-body interactions. These expressions will be 
used in the following sections to analyze the 2D signals obtained by numerical solution of 
equations of motion with finite-envelope pulses. 

We consider a three-pulse sequence where the electric field is given by 

3 

E(r, t) = J2 - tj)^-^) + e . S - {t _ r .) e -^—M] . (i) 

i=y 

Here £f {£j = (£^)*) is the envelope of the positive(negative)-frequency component of 
the j-th pulse centered at Tj, with carrier frequency Uj, polarization unit vector e^, and 
wavevector kj-. The ability to vary the envelopes, polarization directions, durations and 
time intervals, tune the frequencies and even control the phases of optical pulses, provides 
a broad class of possible techniques. By scanning these various parameters it is possible to 
design a multiple-pulse technique^ 1 ^ for a specific application, as is done in NMR.— 
The total Hamiltonian is 

H = H -fi-E(T,t), 

where H Q is the free system Hamiltonian and fi is the dipole operator. The third-order 
contribution to the system's polarization induced by the interaction with the field can be 
written as 

poo /*oo /*oo 

P e4 (r,t)=/ / / dti^tfta^eae^^.^.tO^Cr.t-ts)^^,^-^)^^,^-^-*!), 
Jo Jo Jo 

(2) 

Here, is the third-order response function calculated using time-dependent perturbation 
theory 

4?e 3 e 2ei (M2,tl) = I 3 ([[[£*(*), £es(* - * 3 )] " *3 " h)] , fi ei {t - t 3 - t 2 - h)}) (3) 

where f£ e . (t) are interaction picture operators. is a tensor of rank four, but for clarity we 
hereafter suppress the tensor indices ej. In impulsive experiments where the field envelopes 
are much shorter than the delay periods and the exciton dynamics timescale, the system 
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first interacts with the k x pulse, then with k 2 and finally with k 3 , and the time integrations 
in Eq. (j2J) can be eliminated. We then get 



P^(T t) = R^\t^ t 2 ^ 1 )£^£;^£;^ e i ( ±kl±k2±k 3) r e _i ( ±Wl± ^2±^3)t e «(±^l±^2±W3)*3 e «(±^l±^2)i2 e ±«Wltl 

We shall recast the polarization in the form 

4 

P( 3 \r,t) =Y,Ps(t3,t 2 ,t 1 )e^ r -^ t + c.c, (4) 

s 

where the signal wavevector k s can assume one of the four phase-matching values: ki = 
-ki + k 2 + k 3 , kn = ki - k 2 + k 3 , k m = k x + k 2 - k 3 and ki V = ki + k 2 + k 3 . 

By employing heterodyne detection and controlling the relative phases of the heterodyne 
and incident fields, one can detect the signal field itself (both amplitude and phase) rather 
than merely its intensity (modulus square). In thin samples, propagation effects are negli- 
gible and the signal field is proportional to the polarization field. Hereafter we denote the 
signal field amplitudes generated along k,, as Sj, where j =1, II, III and IV. These are given 
by 

Si (t 3 , h, h) = iRfih, t 2 , ti)£+£+£f e^^^V^^^e - ^' 1 , (5) 
Su (hMM) = ^S ) (t3,t2,ti)^ 2 "^i + e i(wi ^ 2+W3)t3 e i(wi ~ a ' 2)t2 e iu ' ltl , (6) 
Sm {ta, h, h) = iR^{t 3 , t 2 , ^StSte^^-^e^^e^, (7) 
Stv (h, t 2 , h) = iR^ih, t 2 , tl )S+S 2 + S+e l{wi+ ^ + ^ )t3 e l ^ 1+ ^ t2 e lUlt \ (8) 

The exciton level-scheme^ 4 - for the two-band model of semiconductor shown in Fig. l(top 
right ) consists of three manifolds: the ground state (g), single-exciton states (e) and two- 
exciton states (/). The dipole operator only connects the g to e and e to / manifolds. 
Within the rotating wave approximation (RWA), R^ vanishes for this model, the response 
functions Rj [Eqs. (jSJ) to fl2J)] for the other techniques are represented by the double-sided 
Feynman diagrams^ 1 ^ shown in Fig. 1. 
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The sum-over-states expression for is 



#i 3) (^2,ti) (9) 
= z 3 J> a • iJ, ge ,)(e 2 • ^)(e 4 ■ ^ e )(e 3 ■ ^e-C^+r^-r^+^.-r^ 

e,e' 

+ z 3 ^( ei • Mge ,)(e 3 ■ /v 9 )(e 4 ■ ^ e )(e 2 • M Jr^+W^+r>+K-^,l*. 

e,e' 

- z 3 ^ (e a ■ /x ge ,)(e 4 • /x e , / )(e 3 • M/e )(e 2 • / x efl ) e -(-/e'+r /e ')*3-(- ee , + r ee ,) t2+ (- e , ff -r e , s ) il) 

e,e'J 

where u vv > = e v — e v > (u, v' = g, e, e', f ) is the frequency and Y vv i are the dephasing rate of 
the v — > v' transition. The three terms correspond respectively to diagrams (i), (ii) and 
(hi). 

The 2D signal is displayed in the frequency-domain by a Fourier transform of Si (£3, t 2 ,ti) 
with respect to t 3 and t\, holding t 2 fixed 

Sf(Q 3 ,*2,^i) = J J dhdt^hMMY^e^. (10) 

This gives 

<->i C2,"iJ — >_> To 1 1 1 -p irn — I T^F — f 

^ [il 3 - U eg - LOi + L0 2 + UJ 3 + Zl efl J[«l + ^e'g - 0>1 + Zl e /gJ 

| e-^+r ee 0fe( ei ■ ^(e, ■ /y fl )(e 4 ■ /x ge )(e 2 ■ /x eg )g 3 + g+gr 
£^ [fi 3 - u eg - (jOi + u 2 + uj 3 + iT eg ) [pi + cj e / g - + zT e / g ] 
y^ e-^+r ee ,) f2(ei . . . . ^E+E+g- 

^l^f [^3 -w/e' -<^i +o; 2 + o;3 + zT/ e /][Oi +w e ' g - u x + iT e > g ] 

Eq. (ITTi) has various diagonal peaks (fii = fi 3 ) and cross-peaks 7^ f2 3 ). The relative 
contributions of different terms may be controlled by the carrier frequencies, u x , u 2 , and u 3 . 
Spreading the signal in an extra frequency dimension enhances the resolving power of the 
2DCS, compared to ID techniques.- We can further improve the resolution by controlling 
other parameters such as the pulse polarization directions, carrier frequencies and envelopes. 
Other 2D techniques generated in different phase-matching directions and using different 
pairs of time variables (e.g. t 2 and £3 ) provide complementary information^M^ through 
different projections of the response, as will be discussed in Section III. Closed expressions 
for the other 2D signals Sn and Sni are given in Appendix A. 
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III. APPLICATION TO III-V SEMICONDUCTOR QUANTUM WELLS 

A. Single and Two-exciton Resonances 

The dipole selection rules for HH and LH excitons in III-V semiconductor quantum wells 
are shown in Fig. 2. 

For GaAs, the electron effective mass m e = 0.065m (m is free electron mass), while 
the in-plane LH and HH masses are m n LH = 0.206m and m" HH = 0.115m respectively, and 
perpendicular LH and HH masses are mj^ H = 0.094mo and mjj H = 0.34mo respectively. 
The allowed transitions near the bandedge are denoted by right (<x + ) and left (cr~) arrows, 
representing right and left circularly polarized photons, respectively. The corresponding 
transition dipoles ji vc are 



H™ = -^/i a+ = l=no (x+*y) , (12) 

/i " 2C2 = TI^ - = Ti^ ' ^ 

^ = ^ a + = -^o (x+*y) , (14) 

v V4C1 = = (*-*y) , ( 15 ) 

where x and y are unit vectors and fj, is the modulus of the transition dipole. 

We denote each type of exciton by the (hole, electron) spin values involved in the valence 
to conduction band transition. According to Fig. 2, there are two types of HH excitons, 
(— |, — |) and (|, |), and two types of LH excitons, (— |, |), (|, — §)■ The system has 10 
types of two-excitons: three made of two single HH excitons are denoted three made 
of two single LH excitons ( and four made of one LH exciton and one HH exciton 
(f M , mixed). Two electrons (two holes ) belonging to different bands generate a bound 
two-exciton 

[!] ("I."!) + (§>!) = bound/*, 

[2] + : bound / L , 

[3] (-§>-!) + (-!>§) : bound / M , 

W (f,l) + (i-l) : bound / M . 

Other combinations produce unbound two-excitons 
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[5] (-§,-§) + (-f,-|) : unbound f H , 
[6] (14) + (ll) : unbound 
[7] + : unbound f L , 

[8] & -I) + (i -I) : unbound f L . 

The other possible two-excitons (— §, — |) + (|, — |) and (|, |) + (— |, |), are expected 
to form unbound mixed two-exciton /m because the electron motion dominates the internal 
motion of a two-exciton when the holes are much heavier than the electrons. 

In Fig. 3, we show the level-scheme that includes the single-exciton and the four bound 
two-exciton transitions.— 1 ^ The unbound transitions, neglected in Fig. 3 for clarity, will 
be included in the numerical calculations reported in Sec. V. 

The exciton to two-exciton transition dipoles are collectively denoted by /i e ^, where e can 
be either e# (e' H ) or e^e^) and / can be either fi or fu- When all pulses are either 
o~ + or cr~ polarized, the model of Fig. 3 reduces to two coupled two-level systems (dimer) 
whose 2DCS was studied in detail.— 

B. Feynman diagrams for 2DCS 

We first examine the Liouville space pathways^^ for the Si technique. Starting with 
diagram (i) in Fig. 1 we let the states e and e' assume all possible single-exciton states, as 
shown in Fig. 3. This yields the four Feynman diagrams (ia), (ib), (ic) and (id) depicted 
in Fig. 4. 

The corresponding two-dimensional correlation spectrum is shown at the bottom of Fig. 
4. The diagonal peak (ia) comes from HH excitons. (ib) is similarly a diagonal peak for 
LH excitons. The cross peak (ic) describes a pathway with LH excitons during t\ and HH 
excitons during £3. (id) is a second cross peak representing HH excitons during t% and LH 
excitons during £3. Diagram (ii) of Fig. 1 differs from (i) only during the £2 period where 
the system is in the excitonic rather than the ground state. Consequently the corresponding 
peaks (iia), (iib), (iic) and (iid) in the (0 3 , VL{) correlation plot which show the evolution 
during t\ and £3 are the same as (ia), (ib), (ic) and (id), as shown in Fig. 4. 

We next turn to diagram (iii) which involves two-exciton states. Solid (open) symbols 
represent bound (unbound) two-excitons which are red (blue) shifted along ^3 with respect 
to the single-exciton peaks. Running the states e, e' and / over all possible values results in 
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the six diagrams shown Fig. 5. The corresponding 2D spectrum is shown at the bottom of 
Fig. 5. In the total 2D spectrum (Fig. 6(a)) which combines all spectra of Fig. 4 and Fig. 
5, the two-exciton contributions (iiia) to (iiid) appear near the four single-exciton peaks of 
Fig. 4 (ovals and circles). There are four pathways leading to J'm, one to fn and one to /x 
respectively. The two-exciton peaks are either red or blue shifted along fl 3 relative to the 
major diagonal peaks [(en, -eij), (ex,, —ex)], or cross peaks [(e#, —ex), (ex, —en)]- 

To illustrate the merits of 2DCS, we compare the 2D spectra of Fig. 6(a) with a traditional 
spectrally resolved FWM (1DFWM) shown in Fig. 6(b). Using the current notation, the 
1DFWM is given by Si (^3^2,^1), where the two pulse-delays ti and t 2 are held fixed and 
f^3 is varied. The peak denoted HH comes from pathways (ia), (ic) and (iia), (iic). The 
LH peak originates from the other four pathways (ib), (id) and (iib), (iid). In the 2D 
spectra (Fig. 6(a)), (iia) is separated from (iic) and (ia) is separated from (ic) along the Oi 
axis. We find three two-exciton contributions to the 1DFWM in Fig. 6(b) in each side of 
the single-exciton peak HH, and three in each side of the LH peak. These contributions 
are very close and poorly resolved in the ID plot because the typical two-exciton binding 
energies in the III-V quantum wells are 1 to 2 meV. Between e# and ex (3.8 meV for the 
quantum well considered), there are six types of two-excitons congested in the ID plot. In 
contrast, the bound mixed two-excitons (e.g. solid trapezoid and square) and bound HH 
two-exciton (e.g. solid triangle) are well separated along Q x in the 2D plot in Fig. 6(a). 
Moreover, additional separation can be obtained by different projections of the Si signal. 
For example, the overlapping mixed two-exciton contributions such as square and trapezoid 
in Fig. 6(a) (displaced for clarity) can be separated if we plot the 2D spectra Si (fi 3 , Sl 2 , ^i), 
rather than the Si (fi 3 ,t 2) — ^1) in Fig. 6(a) (see details in Appendix B). 

To summarize this section, we have demonstrated how 2DCS can separate different path- 
ways. In Fig. 6(a), coherences among LH and HH excitons correspond to the two major 
cross peaks [(e#, —ex), (ex, — &h)] while the many-body two-exciton correlations correspond 
to the weak peaks around the four major peaks. 2DCS with different choice of variables or 
along different phase- matching directions provide complementary information. 
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IV. HAMILTONIAN AND THE NONLINEAR EXCITON EQUATIONS 

The sum-over-states expressions for 2D signals given in the previous section provide an 
intuitive tool for analyzing 2D techniques. Calculating the signals using these expressions 
requires the single-exciton and two-exciton states and their transition dipoles and relies on 
the impulsive-pulse assumption. This is easily done for Frenkel excitons. An alternative 
way to proceed which is more practical for Wannier excitons is to employ the Nonlinear Exci- 
tonic Equations (NEE)^^^^'^ or the equivalent Dynamics Controlled Truncation (DCT) 
formalism^^^^ to account for the many-body interactions beyond the Hartree-Fock 



\ eve \M&^^Zl&&^ZMMLMZ2M£l To close the infinite hierarchy of dynamic variables, 



the equations of motion are truncated according to the desired order of the laser field: JU&iiUii 
2D signals may be calculated by solving the NEE using nonequilibrium Green's functions and 
the two-exciton scattering matrix— However, in this paper we use direct integration 

of the equations of motion. 

Calculating the 2DCS of semiconductor quantum wells where excitons are spatially con- 
fined in two-dimensions requires an intensive numerical effort. Most computational work 
on the nonlinear optical response of two-exciton correlations focused on HH excitons. Our 
computation time is around 300 times higher than the corresponding ID FWM calcula- 
tion. To make these calculations more tractable, we shall use a multi-band one-dimensional 
tight-binding YaodelM-^^^^^-^^ to describe the excitons and two-excitons of a single 
quantum well. This model can reproduce many spectroscopic observables in quantum-wells 
such as the signs of energy shifts, bleaching and induced absorption and often even their 
relative strengths, and the dependence on the polarization directions of the incident pulses.— 
Comparison of the one-dimensional tight-binding with two-dimensional models can be found 
in Ref.— . Other models have been successfully used as well.— 

We start with the multi-band tight-binding Hamiltonian i s 2±i2&2!LiiLlMiL&& 



describes free band motion, afi^af) are creation (annihilation) Fermi operators of electrons 



H = H K + He + Hi, 



(16) 



where 




(17) 
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in site i from the conduction band c and a^(a^) are creation (annihilation) Fermi operators of 
holes in site i from the valence band v. The diagonal elements Tf? v describe the site energies 
for the electrons (holes) in the conduction (valence) band while the off-diagonal elements 
T5^ represent the couplings between different sites. We adopt the nearest-neighbor tight- 
binding approximation for the electronic coupling, i.e. T^ v = for \i — j\ > 1. The site 
energies of the electrons and holes T^ v are taken to be the mid bandgap energy, E g . 
The Coulomb term He in Eq. (Tl6|) has the monopole-monopole form 

Hc= £ (of + 4 - afV) V* [fh C i ~ > 

ijcvc'v 1 

where the Coulomb interaction is given by 

Vij = U - ^r^. (18) 

\i — ] \ a + a 

This is similar to the Ohno coupling used the Pariser-Parr-Pople Hamiltonian of conjugated 
molecules.— JJ characterizes the interaction strength, do is the spatial cutoff and d is 
the lattice constant. 

Finally, the dipole interaction with the radiation field has the form 

Hj = -E(r,t) P, 

where E (t, r) is given by Eq. (CQ) and P is the interband polarization operator 

p=E[/4^;+c.c.]. (19) 

ijvc 

fjL^j are interband dipoles and = a\a^ are interband coherences. The transition dipole 
matrix elements /i^ c are defined in Eqs. (1121) to (115j) . All optical transitions are diagonal in 
the site indices, i and j. For example, Eq. ( fl2l) reads 

tfP = -^Vo(r + Sij = ^VoSij (x+«y) . 

To make the numerical calculation tractable, we neglect the exciton population dynamics 
and work in the coherent limit where we consider only two types of density matrices involv- 
ing single-exciton and two-exciton respectively.- To third order in the radiation field, the 
equations of motion for the first type of density matrix p?? £ire 2i J 22 1 28 J 5i J 73 J 83 J 84 
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_ i—nV? —v VC 

dt P ^ t ex ^ 

= -Vt c v vc 



jnrm 
n m 



Er .n vc . + V-rj 
x m%Pmj * v i]P'. 



vc 



Fij Fife 



klv'd 



n v ' c '\ n v ' c 'n vc — I r> v ' c \ f> v ' c ' vc 
Plk I Plk Pij [Plk I n lkij 



+ E(f) 



pij 



klv'd 



+e (4 c y{ptypt 



klv'd 



(20) 



where t ex describes exciton dephasing time. B^f? = a£ a?a^ is a two-exciton operator 
whose equation of motion is 



• ryv'c'vc tdv'c'vc 

l d~t lkij u, lkij 

El T31>'dVC I fTTU TDv'c'vC 

\ jm Ikim ' mi Ikmj 

m 

I /tic tdv'c'vc I TH) rji/c'uc\ 
km Imij ' ml mkij J 

+ (V lk + V tj + V lk + V i:i - V u - V kj ) B^ c 
-(yik + Vu-Vu-V^pgptf 

+ (V ik + V lj -V u -V kj )pf k c 'p% 1 (21) 

where £&j is two-exciton dephasing time. 

The total interband polarization is given by Eq. (fl9|) where p"^ is obtained by solving 
Eqs. ( 1201 and fT2~Tj) . To single-out a given FWM signal, we must keep track of the spatial 
Fourier components of the interband polarization contributing to different 2D signals (Sj, 
Sn and Sin). For example, to calculate Si we need the ki component of the interband 
coherences, p™'^. The positive-frequency component (first term of Eq. (JT9l) ) of the total 
interband polarization field in this direction is thus given by 

P [kll fe,t 2 ,ti,t) = w (*8,ta,ti,t). (22) 

ijvc 
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The Si signal (Eq. (llOp ) is finally given by 



Si(^ 3 , t 2 , fli) = 




P [kl] (h, t 2 , h, t)e initl e iQ3t3 dt l( it 3 



ijvc 



(23) 



where p"- depend on t 3 ,t2,h and t through Eq. ( l20i) . as well as pulse envelopes and 
carrier frequencies (Eq. ([1])). are calculated by expanding Eqs. fl20|) and (1211) in the 

various wavevector Fourier components.- We have only calculated the signals for t 2 = 0, 
where the third and second pulses coincide. The order by order expansion of Eqs. fl20|) and 
(f2~Tl) to yield the signal is given in Appendix C. 

To obtain the Si signal, we have solved Eqs. ( 1C3I) . ( 1C4I) to ( 1C6I) for different combinations 
of (£ 3 ,ti). The other signals (Sn and Sni) can be calculated by deriving different sets 
of coupled nonlinear equations from Eqs. (1201) and (pTj) for the relevant spatial Fourier 
components. 

V. NUMERICAL RESULTS 

We have employed the ID tight-binding model to calculate the 2DCS signal from a 
20nm GaAs/Alo.3Gao.7As single quantum well (SQW). 51 For the site energies and carrier 
coupling energies defined in ( ITTj) . we used = 8 meV, T^7-' 2 = 4.75 meV (HH band) and 
T^ 3,4 = 2.52 meV (LH band) to account for the in-plane dispersion of the valence-band 
structure in the quantum well.— ^ The site energies T? = j and T^~ 1,2 are taken to be half 
of the bandgap, E g . We used Gaussian pulse envelopes £^ (t — tj) = exp [— (t — t,) 2 , 
where Sj = 0.6 ps, corresponding to a spectral bandwidth of around 5 meV (FWHM). This 
narrow bandwidth excludes the continuum states and allows to control the relative strength 
of the LH and HH exciton transitions by tuning the carrier frequency, since it is comparable 
to the LH and HH splitting (3.8 meV) of the SQW sample. The pulse power spectra and 
the linear absorption are shown in Fig. 7. Calculational details are presented in Appendix 



Because the oscillator strengths of LH excitons are approximately one third of those of HH 
excitons, we have tuned the carrier frequencies of the three optical pulses uj\ = 0J2 = ^3 = oj c 
to uj c = 3 meV (Fig. 7, dotted), which enhances the LH exciton transitions. The frequency 
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scale in Fig. 7 and all remaining figures are relative to the HH exciton resonance which is 
35 meV below the bandgap. 

To calculate the 2D spectra, we first computed pM for fixed t\ from £3 = to 32 ps 
with 660 time grid points by solving Eqs. ( 1C3I) to ( 1C6I) . These calculations were then 
repeated by varying £1 from 2 ps to 18 ps on a 330-point grid. The computation time 
scales as the number of equations, i.e. 64-iV 4 , where N is the number of sites. We chose 
t ex = 2 ps and t bi = 1 ps, rather than t ex = 4 ps and t bi = 2 ps for the dephasing times of 
excitons and two-excitons as in Ref.— . The numerical simulations converge better for using 
faster dephasing rates but at the expense of the spectral resolution. The current values of 
dephasing times (2ps/lps for excitons and two-excitons) do not show two-exciton peaks but 
merely show shoulders. In all calculations, we employed periodic boundary conditions for 
N = 10. Adding more sites did not have significant effect on our calculations. In Fig. 8, 
we show |Si(f23, £2, t\i )| for ti = 2 ps for different basis size iV =10, 12, 14 and 16. The 
main features of the spectra have converged at iV = 10. Higher energy (> 8 meV) peaks 
represent continuum states. These are not converged at N = 10. However, for the narrow 
pulse bandwidth used (Fig. 7), the contribution of these continuum states is very small. 
Testing the convergence is expensive. A single 1DFWM calculation with 16(18) sites for 
current parameters takes more than 20(32) hours on a single processor (64-bit opteron). 
We thus tested the convergence of 1DFWM signal before launching the 2DCS calculations. 
Fast Fourier Transform of the 660 x 330 grid of P^ kl '(£3, £2, £1, £) with respect to £1 and £3 
gave Si(fi3, £2, — Oi). The typical computational time for one 2D plot for iV = 10 is around 
1000 hours on a single processor (64-bit opteron). 

The calculated |Si(f23, £2, — f^i)| signal with two co-linearly polarized pulses (hereafter 
referred to as XX excitation) is shown in Fig. 9-A (log scale). 

The two diagonal peaks (eg, —en) and (e^, — e^) describe respectively the HH and LH 
excitons. Since the pulse is tuned closer to the LH exciton resonance, the peak (e^, — e^) is 
stronger than (eg, —en)- The two major cross peaks (e#, [pathways (ic) and (iic)] and 
(ei, —gh) [pathways (id) and (iid)] describe the coherences which lead to quantum beating 
among LH and HH excitons. These pathways may be separated using Sj (1^3, 1^2? ^1) &s 
shown in Fig. 13, or the Sn technique (Appendix A). 

To show the two-exciton signatures in Fig. 9-A, we have calculated the 2D spectra 
without the two-exciton contributions (Fig. 9-B) obtained by setting the term Bf k f- Vc in Eq. 
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( 1C3j) to zero. Comparing Fig. 9- A and 9-B, we can clearly see a feature (e) in Fig. 9-A 
for mixed unbound two-excitons. This corresponds to pathways, (iiia') and (iiid). Further 
separation of these two pathways can be accomplished using other techniques as described 
in Fig. 13 or in Appendix A. All the two-exciton peaks appear as shoulders adjacent to 
major single-exciton peaks, or as a broadening of the single-exciton peaks. 

The two-exciton feature (b) in Fig. 9-A can be attributed to the bound LH two-excitons 
discussed in Fig. 6(a) and is visible since the optical pulse enhances the LH excitons as 
compared to HH excitons. Some other two-exciton contributions are also seen in Fig. 9-A. 
For example, a weak shoulder (c), two shoulders (a) and (d) are seen in Fig. 9-A, but not 
Fig. 9-B. As expected, the two contributions (a) and (d) become stronger as the pulses 
are further tuned to the red (see Fig. 10-A). Finally, by comparing the 2D plots (not 
shown) with and without the term including the factor E (t) in Eq. (IC3I) . we found that the 
contribution from Pauli blocking^ is negligible. 

To study the dependence of Si on the pulse polarization configuration, we have calculated 
the spectra with a + a + excitation. These spectra shown in Fig. 9-C are similar to Fig. 9- 
A. However, the signature of bound LH two-exciton, the shoulder (b) in Fig. 9-A, is now 
absent. This is because, according to the selection rules shown in Fig. 3, a + a + excitation 
cannot generate either bound HH two-excitons or bound LH two-excitons. 

Fig. 10-A shows the 2D spectrum when the carrier frequency is tuned to uj c = 2 meV 
(Fig. 7, dashed dot). The 2D spectra calculated without two-exciton contributions is shown 
in Fig. 10-B. 

The effect of tuning is demonstrated by comparing Fig. 10-A (u c = 2 meV) and Fig. 
9-A (uj c = 3 meV). The strongest exciton peak (e^, — e^) in Fig. 9-A becomes the weakest 
in Fig. 10-A, while the weakest exciton peak (e#, —gh) i n Fig- 9-A becomes the strongest 
in Fig. 10-A. Furthermore, because the LH excitons are no longer enhanced for u c = 2 
meV, the signature of the LH bound two-excitons, the shoulder (b) in Fig. 9-A, is absent in 
Fig. 10-A. The two-exciton features (c) and (e) in Fig. 9-A become weaker in Fig. 10-A. 
However, the two-exciton contributions (a) and (d) in Fig. 10-A, become stronger. 

By further red-shifting the pulses carrier frequency to uj c = meV (short dot in Fig. 
7) and using cross-linear (XY) polarization (first pulse is polarized in X direction and the 
second in Y direction), we obtain the 2D spectra in Fig. 11-A. For this excitation condition, 
the pulses are resonant with the HH excitons, and LH excitons are barely excited. Due 
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to the narrow spectral pulse bandwidth, the LH and HH excitons are weakly coupled and 
the (ex,, —en) and (e#, — e^) peaks are very weak. For two-exciton contributions, we can 
clearly see a shoulder (a) for bound HH two-excitons and a broadening (b) for unbound HH 
two-excitons, as compared to Fig. 11-B where two-exciton contributions are excluded. The 
two-exciton contributions are more evident for XY excitation due to strong cancellations 
of the excitonic signal components. 84 Experimental measurements of the XY configuration 
also show the elongation along the Q3 axis, consistent with the calculated spectrum. 

We have further calculated the phase- sensitive 2D spectra rather than their modulus. 
Re(Si(^3, t 2 , — ^1)), calculated with uo c = meV (XX), is shown in Fig. 12-A. The peak 
(en, —&h) resembles that experimentally observed (peak RA of Fig. 2 in Ref.—). The 2D 
spectra in Fig. 12 is plotted in the same way as in Ref.— 

The sum-over-states expression for Si (Eq. pip ) may be used to interpret the phase- 
sensitive spectra in Fig. 12-A. For the two-pulse scheme (t 2 = 0) and colinear excitation, 
we define the real parameter C = e _r "* 2 (ei • n ge ,)(e 2 • /V 9 )( e 4 ' A t ge )( e 3 ' t^eg^t^t^i 1 
where we assume the S 3 £ 2 Si * s re& l- Moreover, for co-linear excitation, the two-exciton 
contributions from the third term in Eq. ( fill) are small and will be neglected. We further 
assume the carrier frequencies lu\ = uj 2 = 003 = uj c . Eq. ffTTj) then gives 
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where u eg = uj eg — uj c and u e i g = u e i g — uj c . We rewrite Eq. (Klj) as 

S (3) (Q t Q x y> 2C [(^3 - U eg ) - jV eg \ [gl + Cbefg) - iT elg \ 

' 2 ' £ [(^-^ 9 ) 2 + ry {{n 1+ Q e ,gf + Yi lg ] 



and thus 
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s\ 3 \n 3 ,t 2 ,-n 1 )]=Y ^^-^ (^ + ^ 9 )-2CT eg r e/g 



For the excitation of only HH exciton as in Fig. 12-A, Eq. (1261) is simplified as 
Re 



Sf\(l 3 , fa, -n,)l = r 2C 7 + °™> ~ f"" 1 , , (27) 
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where lohh is the HH exciton transition frequency relative to the carrier frequency and Thh 
is the exciton dephasing rate. 

Eq. (1271) accounts for the main features of Fig. 12-A. The denominator gives the reso- 
nance peak at Q3 = ujhh and Q\ = —Ohh, while the numerator gives the positive/negative 
feature. When neglecting Tj IH , we see that if f2 3 goes from below to above u>hh, the signal 
will change sign for fixed Similarly, when Qi is scanned across -ljhh, the signal will 
also change sign for fixed Q3. In the region where the signal changes signs, steep gradient 
occurs in the form of the dense contour lines in the center of the blue/red peak shown in 
Fig. 12-A. The zero line along which the steepest local gradient occurs must pass through 
the peak center (uhh, —&hh) as at this point Eq. ( l27j) vanishes when T 2 HH is neglected. 
However, for finite T 2 HH , Eq. ( 1271) is finite at (ujhh, —&hh) and thus the zero contour line 
with steepest local gradient, or the line passing through the center of the negative/positive 
(blue/red) peak will not pass through the point (ujhh, —&hh), which is on the diagonal line. 
Thus, due to the finite dephasing, the peak center will be off the diagonal line. Because 
T 2 HH is dominated by Coulomb interactions beyond first-order (carrier-carrier scattering), 
the asymmetry of the peaks and their offset from the diagonal lines may provide insights on 
these many-body correlations. The above analysis applies only to the excitation of a single 
level, e.g. HH excitons. The picture will become more complicated if a multi-level system 
is excited. 

The calculated 2D spectra (real part) with carrier frequency u c = 2 meV (XX) shown in 
Fig. 12-B share some similarities with recent experiments^, especially for the main peak 
(en, — e#). However, it is difficult to produce 2D spectra with all four peaks observed in the 
experiment. In our calculations, we use a different SQW as compared to the experiment. 
A two-pulse scheme is employed in our simulations while the experiment was done with 
three pulses. In addition, the 2D spectra are very sensitive to the carrier frequency of 
the narrow-spectra optical pulses and to the parameters such as site energies, hopping and 
dephasing rates for the HH and LH excitons. For example, by tuning the pulse frequency to 
the blue (u c = 3 meV (XX)), the real signal (Fig. 12-C) changes considerably compared to 
Fig. 12-B. The experimental data also show strong dependences on tuning of the excitation 
pulses as well as their strength.—^ Finally, incorporating the population dynamics and 
continuum states will be required to fully reproduce the experimental spectra. 
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VI. CONCLUSIONS 



As in NMR, optical 2D techniques supplement traditional ID techniques by revealing mi- 
croscopic couplings, coherences and many-body correlations. Using the Liouville pathway 
analysis, we have predicted the exciton and two-exciton contributions to 2DCS in a single 
quantum well with LH and HH exciton couplings. Numerical simulations reproduce the ma- 
jor diagonal peaks and cross-peaks in a 2D spectra. We also obtained several features from 
two-exciton contributions. Different 2D techniques provide complementary information. 

Some features predicted by the sum-over-states expressions, for example the bound mixed 
two-excitons, are not seen in our simulations. This may be due to the small binding energies 
of these bound two excitons and to the excitation conditions used. Future extension of the 
current two-pulse scheme to three-pulses^ and exploring the independent 2D spectra will 
be of interest. 
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GM59230-05) is gratefully acknowledged. We wish to thank Drs. T. Meier, D. Abramavi- 
cius, M. M. Dignam for most useful discussions. 



APPENDIX A: THE S n AND S m TECHNIQUES 

The response function for the kn = ki — k2 + k3 signal is given by 



flfiW2,*i) (Al) 
= * 3 5> 4 • Av)(e 3 • AV s )(e 2 • M ffe )(e 1 • ^e-^.+^-r^-^+r^ 

e,e' 

+ i 3 J> 2 • AvXeg • /v 9 )(e 4 • Mge )( ei • » eg )e-^ +r es )t 3 -(^+r ee ,)t 2 -(^ + r es ) tl 

e,e' 

e,e' / 

where the three terms are represented by the Feynman diagrams (iv), (v) and (vi) of Fig. 
1 respectively. Substituting Eq. I I All) into Eq. (jBJ) gives time-domain signal Snfa, t^, ii). 
The 2D correlation plot is then obtained by Fourier transforming Sn(ts, t\) with respect 
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to ti and t 3 , holding t 2 fixed, i.e., 

S^(Q 3 ,hM = / f dtzdtxSnfahitJe^e* 1 ^. (A2) 

We get 

S®(n 3 ,t 2 M (A3) 
e-r 99 *2(e 4 • Av)( e 3 • /V 9 )( e 2 • M 9e )( e i ' /O^^i" 



e,e' 
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EE 
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, ,. [fi 3 - w/ e ' + ^ 3 - ^ 2 + + ir /e /] [fix - u; e3 + W\ + «r eg ] 
The response function for the k m = k x + k 2 — k 3 signal is given by 
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where the two terms correspond to the Feynman diagrams (vii) and (viii) in Fig. 1. Sub- 
stituting Eq. (IA4B into Eq. ([7]) gives time-domain signal Sni(t 3 , t 2 , ti). The 2D correlation 
spectra are obtained by Fourier transform Sm{t 3 ,t 2 ,t 1 ) with respect to t 2 and £ 3 , holding 
ti fixed, 



This gives 
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Sn and Sin techniques provide complementary information to the Si technique. Sin 
technique is particular sensitive to exciton-exciton coupling. The 2D spectra of techniques 
Sn and Sin can be similarly analyzed as in the Si technique. 

APPENDIX B: ALTERNATIVE 2D PROJECTION OF THE Si SIGNAL 

The bound mixed two-excitons from different pathways (e.g. solid trapezoid and square) 
not resolved in Fig. 6(a) may be separated by displaying the Si 2D spectra with different 
variables such as (f^, Q2, ti) rather than (f2 3 , t 2 , — ^1), or by employing other 2D techniques 
(Appendix A). In Fig. 13 we display Si(fi 3 , fi 2 , h). Panels (a) and (b) correspond to 
the Feynman diagrams in Fig. 4 and Panel (c) corresponds to the Feynman diagrams in 
Fig. 5. Panel (d) shows the total spectrum from the diagrams of Figs. 4 and 5, where 
the single-exciton contributions along ^3 (from Fig. 13(a) and (b)) are represented by two 
open vertical ovals. 81(^3,^2,^1) can resolve peaks indistinguishable in Si (Q3, — 
For example, (ic) and (iic) (both denoted by a horizontal oval) overlap at (en, —zl) in Fig. 
6(a). Similarly, (id) and (iid) (both denoted by a vertical oval) overlap at (e L) —e H ). The 
two-exciton contributions denoted by a solid trapezoid and a solid square in Fig. 6(a) also 
overlap. The same holds for the two-exciton contributions denoted by a solid hexagon and 
a solid inverted triangle in Fig. 6(a). However, all of these overlapping couplings can be 
easily resolved by 81(^3,^2,^1). In Fig. 13(a), pathways (ia) and (ic) both contribute to 
the peak (Q 3 ,Q 2 ) = (e#,0) (displaced in the figure for clarity). However pathway (iic) 
now gives a peak at (en, en — e£) (Fig. 13(b)) and is thus well separated from (ic), in 
contrast to Si (fi 3 , t 2 , — f2i) (Fig. 4) where these two pathways are overlapped. Similarly, 
pathway (iid) appears at (e^,eL — e#) (Fig. 13(b)) and is thus separated from (id), which 
appears at (ei,0) in Fig. 13(a). Thus Si (^3, ^2, ii) (Fig. 13) provides a complementary 
information by separating different pathways which overlap in Si (f2 3 ,t 2 , — ^1) (Fig. 6(a)). 
All two-exciton contributions to Fig. 5 are plotted in Fig. 13(c). Fig. 13(c) also provides 
complementary information to Si (0 3 , t 2 , — ^1) (Fig. 6(a)). For example, the mixed two- 
excitons from two pathways (solid square and a trapezoid) overlap in Fig. 6(a) but are well 
separated in Fig. 13(c). The same holds for the solid hexagon and solid inverted triangle. 
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APPENDIX C: SELECTING SPATIAL FOURIER COMPONENTS OF 
INTERBAND POLARIZATIONS 

In this appendix, we present the equations of motion used for calculating FWM signals 
for a two-pulse scheme (t 2 = 0). According to Eq. ( 1201) . the kj spatial component of 
interband density matrix, p" c '^, satisfies 
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dt J 
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To first order in the optical field, the single exciton density matrix p contains either ki or 
k2. Invoking the rotating wave approximation (RWA) and by a perturbative expansion^* 
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to 3rd-order in the optical field, Eq. (1C2j) becomes 
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Apart from the density matrix element p v ^ kl+k2+k2 ] 5 gq (]C3|) also contains the matrix 
elements P ; c[kl] , P ; c[k2] and <^ c:[k2+ka] . Thus to solve for p; c:[ - kl+k2+k2] (t 3 , t 2 , t u t) at time 
t, we should calculate p^ kl ' (t), pjf k2 ^ (t) and B^cvcfa+ka] (t) simultaneously. Following 
the similar procedure of deriving Eq. (1C3j) . we obtain the following equations for solving 
ftf l] (*), pf^ (f) and <^ c:[k2+k2] (t): 
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Solving the set of coupled equations (1C3I) . (1C4[) to (1C6j) gives p-^ ki+k2+k2 '(t 3 , t 2 , ii, t). 
Substituting p^ c '' _kl+k2+k2 '(t3, t 2 , ti, t) into Eq. (j2"51) results in the Si 2D signal. Equations 
for the Sn and Sni 2D signals can be derived in a similar way. 

APPENDIX D: NUMERICAL CALCULATION OF LINEAR ABSORPTION 

From Eqs. ( Tl9|) and ( l22l) . the total polarization along the ki direction is written as 



P(t)=pW(t3,t 2 ,tl,t)+C.C. 

= E^ ;[kll ( t 3^2,t 1 ,t)+C.C. 



-iiOnt 



= p(t)e- lU3at + c.c, (Dl) 

where p (t) is defined as the slowly- varying portion relative to the bandgap frequency, u g . 
Performing Fourier transform of the polarization relative to u g gives 
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Similarly, performing Fourier transform of the j—th optical pulse (Eq. (pQ)) relative to 
the bandgap frequency u g gives (spatial portion is omitted) 
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where E opt (w) is approximately the Fourier transform of pulse envelope £j~ because we 
consider only resonant excitation where Uj is close to uj g . Finally, the linear absorption is 
calculated as 36,94 
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where n' is the average, frequency-independent refractive index of the quantum well, eo is 
the vacuum permittivity and c is the speed of light. 
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Figure Captions 

Fig. 1. Feynman diagrams for the possible coherent 2D spectroscopic techniques for the 
exciton level scheme given in the top right. \g) is the ground state, |e) is the single exciton 
manifold and | /) is the two-exciton manifold. /i 9e and /z e / are the corresponding transition 
dipoles. [From Ref.-] 

Fig. 2. Selection rules for III-V semiconductors. c\ and C2 are two conduction bands 
whose spins are s z = — \ and s z = | respectively. v\ and t>2 are HH valence bands with 
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spins s z = — | and s z = | respectively. v 3 and t> 4 are LH valence bands with spins s z = — | 
and s z = | respectively. 

Fig. 3. Single-exciton and bound two-exciton optical transitions in III-V semiconductors. 

Fig. 4. Top: the Feynman diagrams derived from Fig. l-(i) and (ii). Bottom: the 
schematic 2D spectrum. 

Fig. 5. Top: the Feynman diagrams derived from Fig. l-(iii). Bottom: the schematic 
2D spectrum. 

Fig. 6. (a) The total schematic 2D spectrum of ki = k] + k 2 + k 3 (diagrams (i), (ii) 
and (hi)), (b) The 1DFWM signal |Si (f2 3 ,t 2 ,ti)| displayed versus f2 3 for fixed t 2 and t\. 

Fig. 7. (solid): Simulated linear absorption of the single quantum well (see text for 
parameters). Also shown are optical pulse power spectra |E opt (o;)| 2 (Eqs. ( 1D3I) ) for u c = 3 
meV(dotted), u c = 2 meV(dash dot) and u c = meV (short dot). 

Fig. 8. Convergence of the modulus |Si(0 3 , t 2 , ti)\ with basis size. The 0-6 meV region 
has converged for iV = 10. The higher energy continuum has not converged. 

Fig. 9. (Color online) (A) The modulus |Si(0 3 , t 2 , — ^i)| for co-linearly polarized pulses 
(XX) (uj c = 3) meV; (B) same as (A) but excluding two-exciton contributions; (C) same as 
(A) but with a + a + excitation. 

Fig. 10. (Color online) (A) The modulus |Si(r2 3 , t 2 , — f^i)| calculated with co-linearly 
polarized pulse excitations (u c = 2 meV); (B) same as (A) but excluding the two-exciton 
contributions. 

Fig. 11. (Color online) (A) The modulus |Si(f2 3 , i 2 , — ^i)| calculated with cross-linearly 
(XY) polarized pulse excitations (u c = meV); (B) same as (A) but excluding the two- 
exciton contributions. 

Fig. 12. (Color online) Re(Si(0 3 , — fii)) calculated with XX excitation for (A) u c = 
meV; (B) uj c = 2 meV; and (C) uj c = 3 meV. 

Fig. 13. The Si (Q 3 , f2 2 , ii) signal, (a) Diagrams (ia), (ib), (ic) and (id), (b) Diagrams 
(iia), (iib) and (he), (c) Diagrams (iiia) to (hid), (d) Total spectrum. 



(0 | 8 )(g 




(iv) | g )(g 







k 3 


V) 




(e'\ 


k. » 

1 






(v) 


„ \/ „ 

g ){g 
























s)<g 




(vi) 1 1 


/ \ 




kT\ 






1/) 






kf^ 








(vii) 



ki l*X*l 

k n =kj-k 2 +k 3 



I/) 

k ^\g >0 

(viii) 

k^ 

\f) 

k i | g )(g 
k m -k, -k 2 +k 3 



Fig. 2 Yang 




Fig. 3 Yang 




* H )( g H 

\f H 

>0 

(iiia) A A 



l*><* 
(iiito fu 



I/; 



a; 



-1, 

*><* 

(iiia) # O 

k 4 kx^i 

|/m) 



i*><« 

(iiic) ■ □ 



\ e L)( e L 



>0 

(iiib) 4 o 

K )( e H | 

|/m) 



l*><* 
(iiid) ▼ V 











WW \JJ 


♦ o 




A A 


#▼ vo 










Fig. 6 Yang 



Fig. 7 Yang 

i ■ 1 ■ 1 ■ 1 ■ 1 n 




Photon Energy (meV) 



Fig. 8 Yang 

m — ■ — i — ■ — i — ■ — i — ■ — i — ■ — i — ■ — i — ■ — i — ■ — i — ■ — n 




l_i . i . i . i . i . i . i . I i_l . l_J 



-4 -2 2 4 6 8 10 12 14 

^ 3 (meV) 



a si 

*l 1 t- a *i q 
t- n » t- «r t- 




(A9LU) [ U- 



LL 





Fig. B-1 Yang 



